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Electric field gradient effects induced by an asymmetrically in-plane shaped GaAs/AlGaAs 
modulation-doped structures of various design are investigated within 4-300 K temperature range, 
ft is demonstrated that current -voltage characteristics of such structures at low, 4-80 K, tempera- 
tures exhibit well-pronounced asymmetry arising due to a presence of two different gradients of the 
electric field in a two dimensional electron gas. This phenomenon is caused by both, different accu- 
mulation of two-dimensional electrons due to asymmetrical shape of the structure and nonlocality 
in the electron drift velocity. Experiments are illustrated by a phenomenological model and Monte 
Carlo simulation. Possible applications of the effect to detect electromagnetic radiation of GHz and 
THz frequencies are discussed as well. 
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I. INTRODUCTION 

As it is known, in modulation-doped (or selectively- 
doped) structures free electrons, spatially separated 
from ionized donors, are confined to a narrow potential 
well forming the so-called two-dimensional electron gas 
(2DEG)i'^. Along the z axis, i. e. perpendicular to the 
2DEG plane, the confinement potential has a triangu- 
lar shape, and the spectrum of the 2DEG is composed 
of a sequence of quantum subbands.'^ The shape of the 
confining potential has influence only to the carrier tran- 
sitions between different subbands; meanwhile for exci- 
tations within each subband - so-called intraband exci- 
tations - this effect is very small. Perpendicular to the 
z axis, i. e. in the 2DEG, or x — y plane, the carriers 
display very high electron mobilities exceeding values of 
2x10^ cm^/(V-s) at liquid nitrogen and 1x10^ cm^/(V-s) 
at liquid helium temperatures.^ This fundamental prop- 
erty of these structures allowed to invent HEMTs - high 
electron mobility transistorsi^ 

The carrier transport in 2DEG layers at low electric 
fields is described by the Ohm's law together with the 
Einstein relation for the diffusion coefficient, i.e. 



v^ fi[E - (fcr/e)(l/n)Vn] 



(1) 



where v is drift velocity, fi stands for electron mo- 
bility, E is the electric field and n is the free carrier 
concentration; e is the electron charge, k and T de- 
note Boltzmann constant and lattice temperature, re- 
spectively. This so-called "drift-diffusion equation" de- 



scribes the quasi-classical situation when the quantiza- 
tion of the electronic states-induced effects are not im- 
portant. On the other hand, it does not take into account 
the hot-electrons phenomena when the applied in-plane 
electric field is strong enough to heat the electrons up to 
the energy exceeding their equilibrium valuci^ 

The electric field in modelling, as a rule, is assumed 
to be uniform. Although this assumption is very con- 
vienent for calculations, however, in many cases, in or- 
der to obtain the real picture of physical processes, it is 
necessary to include spatial variation of the carrier con- 
centration or/and the electric field along the channel. In 
particular, when it becomes comparable with the range 
of a mean free path. Under these circumtances, effects 
such as velocity overshoot begin to predominate defining 
thus the performance of the device. Such an approach 
is very important, for instance, in understanding the be- 
havior of devices containing asymmetric channel profile 
produced by a relevant doping. More specifically, in a 
quarter-micron n-type silicon metal-oxide semiconductor 
field effect transistors (MOSFET) with asymmetric chan- 
nel profile, formed by the tilt-angle ion-implantation af- 
ter gate electron formation, allows one to achieve high 
current drivability and hot-electron reliability!^ For ex- 
ample, in 0.1 /xm gate length asymmetric n-MOSFET 
structures this technological innovation allows to attain 
higher electron velocity in comparison with conventional 
devices of such type 4 

In this article, we report on experimental and the- 
oretical investigation of the electric field gradient- 



induced effects due to asymmetrically in-plane shaped 
GaAs/AlGaAs modulation-doped structures. We show 
that current-voltage characteristics of such structures at 
low, 4-80 K, temperatures exhibit pronounced asym- 
metry. The physics behind is attributed to a two- 
dimensional bigradient effect which is equivalent to the 
phenomenon observed earlier in a bulk asymmetricaly- 
shaped semiconductor s^'^" . We demonstrate that de- 
pending on the values of the in-plane electric fields and 
their gradients, the effect can reveal itself as a result 
of different distribution of accumulating two-dimensional 
electrons due to the asymmetrical shape of the structure, 
and/or the exhibition of nonlocal drift velocity which be- 
comes pronounced in a different manner due to the pres- 
ence of two different gradients of the in-plane electric 
field. 

The paper is organized as follows. In Sec. II we present 
the design of the GaAs/ AlGaAs modulation-doped struc- 
tures, their electrical parameters and geometry features; 
we also describe briefly used measurement techniques. 
Section III reports on experimental results obtained in 
GaAs/AlGaAs modulation-doped structures of various 
designs and at different lattice temperatures. Section IV 
is devoted to theoretical models and illustrates the con- 
cept of the electric field gradients-induced phenomenon 
- the bigradient effect with special emphasis on mani- 
festation of the electron drift velocity. In Sec. V pos- 
sible applications of the effect for the sensing of elec- 
tromagnetic radiation within GHz-THz frequencies are 
disccused, while in Sec. VI features of three-dimensional 
vs. two-dimentional effect are compared. Finally, con- 
clusions are given in Sec. VII. 



II. SAMPLES AND MEASUREMENT 
TECHNIQUES 

Two types of modulation-doped structures of 
GaAs/Alo.25Gao.75As (structure 2DEG-A) and 
GaAs/Alo.aGao.TAs (structure 2DEG-B) were grown 
by molecular beam epitaxy technique. Their design 
parameters are given in the caption of Fig. 1. 

Electron sheet density, n2D, and low- field mobility, /x, 
at different, 300 K, 77 K and 4.2 K, temperatures for 
structure 2DEG-A are the following: 5.5x10^^ cm^^ and 
4700 cmV(V-s), 1.9x10" cm'^ and 190600 cmV(V-s), 
1.9x10^^ cm~^ and 2x10^ cm^/(V-s), respectively; 
while for structure 2DEG-B: 5.6x10^^ cm~^ and 
8000 cmV(V-s), 4.3x10" cm'^ and 116000 cmV(V-s), 
4.3x10^^ cm~^ and 1x10^ cm^/(V-s), respectively. 

The wafers were then processed into asymmetrically- 
shaped samples of different geometry which was changed 
varying the width w of the neck. We have used the optical 
lithography for the samples with w — 12 /im. The meza 
with height, d, of 2 /xm in these samples was fabricated 
by wet etching. In order to fabricate samples with a 
narrower neck, respectively 5 /im, 2 /im, and 1 /im width, 
we have employed electron-beam lithography and shallow 




FIG. 1: The shape of the studied structures placed on semi- 
insulating substrate. Characteristic dimensions of the struc- 
ture are the following: L=50Q fim; D=250 ^m, A=50 ^m, 
W—100 nm. White color denotes active part containing 
2DEG which is shown schematically as a black sheet; grey col- 
ored parts depict Ohmic contact areas of length C=100 /im. 
Layer sequence of GaAs/ AlGaAs modulation-doped struc- 
tures (from the top): structure 2DEG-A and structure 2DEG- 
B, respectively: 20 nm i-GaAs cap layer; 80 nm Si-doped, 
layer of Alo.25Gao.75As and 60 nm Si-doped, 
of Alo.3Gao.7As; undoped spacers, 45 nm 
Alo.25Gao.75As and 10 nm Alo.3Gao.7As; 1000 nm and 600 nm 
of i-GaAs; smoothing superlattice - twenty and six periods of 
9 nm AlGaAs/1.5 nm-GaAs layers; 0.5 /im and 0.6 nm layer 
of i-GaAs; semi-insulating substrate. 



wet etching; the meza height in this case was 300 nm. 
The Ohmic contacts were produced by a rapid annealing 
procedure of evaporated Au/Ge/Ni compound. 

The wired samples were then mounted into a two-stage 
closed-cycle helium cryostat, and I — V characteristics 
were measured within 4-300 K temperature range. In 
microwave experiments performed at 10 GHz frequency, 
we have used magnetrons delivering pulses of 1.5-5 ^s 
duration with the repetition rate of 35-40 Hz; the sam- 
ples were then placed into rectangular waveguides, and 
the microwave radiation-induced voltage arising over the 
ends of the sample was measured by an oscilloscope; the 
experiments were performed at room and liquid nitro- 
gen temperatures. In the terahertz (THz) range, 0.584- 
2.52 THz, the source of THz radiation was an optically- 
pumped molecular laser operating in the continuous wave 
regime; the illumination into the sample, located in the 
closed-cycle helium cryostat behind the THz transparent 
high-pressure polythene window, was focused by spher- 
ical mirrors. As in the microwave case, the external 
radiation-induced voltage arising over the ends of the 
sample was recorded by lock-in amplifier at a chopper 
frequency of 187 Hz. 



III. EXPERIMENTAL RESULTS 

Figure 2 shows the experimental results oi I —V char- 
acteristics measured in structure 2DEG-A in the sam- 
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FIG. 2: (Color online) I — V characteristics of the asym- 
metrically shaped samples with different neck width at liq- 
uid helium temperatures. The inset depicts the measurement 
scheme: the contact B was always grounded, positive (filled 
symbols) and negative (empty symbols) voltage was applied 
to the contact A. 



FIG. 3: (Color online) I — V characteristics of the asymmet- 
rically shaped sample vs. temperature. The neck width is 
2 ^m. Inset shows the change of asymmetry varying the lat- 
tice temperature. Measurement scheme is the same as in the 
previous plot. 



pies with different neck widths recorded at liquid helium 
temperature. (The experimental data of the structure 
2DEG-B are equivalent and therefore are not presented 
here) . 

As one can see, the salient feature of all the measured 
I — V characteristics is an asymmetry becoming well- 
procounced above 0.1 V. Moreover, the character of the 
characteristics strongly distingtive to the reduction of the 
neck width: For the sample with neck width of 12 /im two 
parts in the curve can be distinguished - the initial rising 
part (below 0.3 V) which bends with voltage and above 
0.7 V transits to the second, saturated one. With the 
decrease of the width w the saturated part transforms to 
a boosting region with highly expressed asymmetry. A 
complex behaviour is observed in the 5 /Ltm width sample: 
careful look can resolve initial rising part below 0.15 V, 
rudiments of saturated part within 0.15-0.6 V and al- 
ready emerging boosting region above 0.6 V with spread 
influence on all the curve. 

It should be pointed out that the asymmetry is 
strongly sensitive to the temperature variation. This is 
illustrated by the experimental data for 2 /im neck sam- 
ple given in Fig. 3. It is seen that at low temperatures the 
asymmetry is rather large - if the coefficient of asymme- 
try of the I —V characteristics is defined as (// — Ir)/If 
- at liquid hefium temperature under 0.8 V it is 0.506, 
at 80 K its value is about 0.369, meanwhile at room tem- 
perature the asymmetry is nearly invisible. To resolve 
the effect, the experimental data should be multiplied by 
a factor of one hundred. 



IV. THEORETICAL APPROACH AND 
PHYSICS BEHIND 



Electric-field gradient-induced effects: carrier 
accumulation 



To understand the physics behind the observed phe- 
nomenon, we have followed the same ideology as in the 
bulk case of asymmetrically-shaped structures: We ex- 
plore an asymmetrically shaped two-dimensional struc- 
ture by a solving a coupled system of one-dimensional 
Poisson and current-flow equations neglecting quantum 
effectsi'i^: 



-^±[S{x)E{x)]^^^-^{n{x)~N,), (2) 
o {X ) dx 6 

/ = {\e\n{xUE)E{x) + \e\D{E{x))^[nix)]}S{x), (3) 

ax 

where x denotes the position along the symmetry axis 
of the asymmetrically-shaped sample, S{x) is the cross- 
section area of the sample, n{x) is the electron con- 
centration, Nfi is the donor concentration (used here 
as a free fit parameter), E{x) denotes the electric field 
strength, e is the dielectric constant of the sample, D{E) 
denotes the diffusion coefficient, and fJ,{E) is the elec- 
tron mobility. The diffusion coefficient may be set to 
a constant value as its dependence on the electric field 
has no appreciable influence on the results. As it is 
known, mobility fJ,{E) changes with the electric fieldi^i 
fi{E) = ^QEcino)/[Eix) + Ec{^J.o)], where Ec = Wc/mo, 
and Vc is the cut-off velocity equal to 10^ cm/s. 

For a given value of the current /, one can employ an 
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tude smaller than that at 4.2 K temperature in the given 
voltage range. The temperature effect is also confirmed 
by theoretical calculations (data are not given here). 

Therefore, the presence of two different electric field 
gradients in the hot two dimensional electron regime, al- 
lows one to explain the asymmetry oi I ~V characteris- 
tics at low temperatures. The phenomenon can be called 
thus the two-dimensional bigradient effect. 

It is reasonable that the reduction of the neck should 
cause stronger asymmetry in I — V curves. Experimen- 
tally it is clearly illustrated in Fig. 2. However, if one 
compares the behavior of the curves, for instance, for the 
5 /im neck sample with a corresponding theoretical es- 
timates (Fig. 5), evident differences can be indentified: 



FIG. 4; (Color online) Calculated electric field and carrier 
density in the sample with 5 /^m neck width under / — 
0.03 mA at liquid helium temperature. 



iterative scheme to solve self-consistently the system and 
calculate the voltage. 

Results of the electric field and relative carrier density 
calculations are given in Fig. 4. The model calculations 
do give two different gradients of the electric field induced 
by the asymmetrical shape of the structure depicted in 
Fig. 1. One must note further that the distribution de- 
pends on the current direction, and this difference is ap- 
preciable only within several microns from the neck: The 
electric field peaks are sharper in the case of the reverse 
current, when the A-contact of the sample is under a 
positive bias. 

Hence, the in-plane geometrical asymmetry induces 
two different gradients of the electric field in the sam- 
ple. As a consequence, different spatial accumulation of 
two-dimensional electrons in the vicinity of the neck can 
be a reason for asymmetric I — V curves. This prediction 
is confirmed by the theoretically calculated I — V char- 
acteristics and their comparison with the experimental 
data presented in Fig. 5. As one can see, the measured 
and calculated I — V curves in 12 /xm neck width samples 
behave similarly: in low electric fields, where Ohm's law 
is still valid, the current values are the same and do not 
depend on the polarity of the applied voltage. However, 
with the increase of voltage - when sample is tuned to 
a strong field regime - the electrons become hot, they 
accumulate differently in the vicinity of the neck due 
the asymmetrical shape of the structure. As a result, 
the asymmetry in the I — V characteristics is observed. 
Since this phenomenon relies on hot-carriers, it is reason- 
able that is temperature-sensitive: the decrease of lattice 
temperature will make the effect much more pronounced 
since the carrier heating at low temperatures is stronger. 
This is clearly seen in Fig. 3: the difference between the 
forward and the reverse-bias curves, for instance, at room 
temperatureJi is approximately by three orders of magni- 
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FIG. 5: Comparison of the measured (top plot) and calculated 
1 — V (bottom plot) characteristics in the asymmetrically- 
shaped samples with different - 12 ^m and 5 /im - neck widths 
at liquid helium temperature. 



the theory does not follow the experiment anymore nei- 
ther in curve shape nor in asymmetry value. Therefore, 
in order to understand the origin of the behavior, a more 
complete theoretical model must be applied. 



B. Electric-fleld gradient-induced effects: Influence 
of nonlocality in drift velocity 

Althought the equation (1) is the backbone of elec- 
tronic device simulations, in the current situation one 
needs to be extended taking into account the electric 
field gradient-induced term,- The electron average veloc- 
ity can be then expressed using the extended/ augmented 



drift-difFusion modeli^ 



S dE Ddn 
E dx n dx 



(4) 



where Vh is the drift velocity in a homogeneous electric 
field, D denotes diffusion coefficient, S depics the length 
coefficient. The latter parameter is in the order of elec- 
tron mean free path and determines the severity of the 
velocity overshoot effect which is pronounced when the 
variation range of the electric field is equal to or shorter 
than S. Also, it is worth noting that in the given equation 
we have assumed that nonlocality in diffusion coefficient 
is negligible and has no observable influence on the effect. 

A point of departure for analyzis is an estimate of a 
mean free path (A) of two dimensional electrons at liquid 
helium temperatures. As electron mobility in the studied 
structures is of about 1x10^ cm^/(V-s), this results to the 
value of 



A 



e 



1 fj,m, 



(5) 



where 77i*=0.068me, the latter is the free electron 
mass. Meanings of other symbols are the same as in 
previous formulae. 

The second step in considering the data is an estimate 
of the distribution of the electric field in the vicinity of 
the neck. Voltage value of 0.25 V was chosen as a refer- 
ence point, since experimentally observed asymmetry of 
the I — V characteristic starts to be strongly expressed. 
Theoretical calculations using the set of equations (2 - 
3) are given in Fig. 6. As one can see, the shape of the 



asymmetry becomes more pronounced. Futhermore, dis- 
tribution of the electric field is squeezed to the narrower 
range. If for 12 /^m width sample the electric field ex- 
tends over 1.7 /ini, in the 5 /im and 2 /xm necks it is 
of about 1.2 /im and 0.9 /im, respectively. One can as- 
sume therefore that in 5 /im neck samples the nonlocality 
can start to predominate since the change of the electric 
field along the structure is comparable with the electron 
balistic length. 

Calculations illustrating the effect of nonlocality in 
electron drift velocity on I — V characteristics are given 
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FIG. 6: (Color online) The dependence of the calculated elec- 
tric field distribution in the vicinity of the asymmetrically- 
shaped sample's neck for two opposite currents at 4.2 K. The 
background is shaded in the form of the sample. Note the 
drastically increased asymmetry in the electric field distribu- 
tion with the reduction of the neck width. 

electric field distribution strongly varies with the re- 
duction of the neck width: it becomes sharper and the 



FIG. 7: Top plot: Calculated 1 ~ V-characteristics of 5 fim 
neck sample using different length coefficient S. Note the in- 
crease of the asymmetry with 5 in the dependences. Bottom 
plot: Experimental results in 5 /im neck sample given for 
comparison. 

in Fig. 7. One can see that increase of S has obvious ef- 
fect in 5 /im neck sample - even small increase in length 
coefficient (in the range of 0.1-0.2 fim) stimulates the 
raise in the asymmetry. It is evident that calculations 
for the 5 /im are very similar to the given experimental 
results. It is deserve to remark that in 12 fim. neck sam- 
ples such the infiuence is not observed (these theoretical 
data are omitted here). 

The inherent feature of the nonlocal effects is a mani- 
festation of the velocity overshoot. To evidence the lat- 
ter we have performed the Monte Carlo calculations of 
the drift velocity in a "frozen" electric field taken from 
a solution of the coupled set of Poisson and current-flow 
equations (Eq. 2-3) as described previously. We have 
simulated the drift of Monte Carlo particles along the 
x-direction in the given electric fleld distribution assum- 
ing that the particles enter the modelling region hav- 
ing the Boltzmann distribution at the lattice tempera- 
ture. The standard three-valley GaAs model is applied 
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FIG. 8: (Color online) Monte Carlo calculations of the elec- 
tron drift velocity in "frozen" electric field in 5 fim neck sam- 
ple. Note oscilations in the drift velocity caused by optical 
phonon emission. The vertical dashed line indicates the neck 
position. Background is decorated schematically as the geo- 
metrical shape of the structure. 



in the simulation-'^^ with an impurity scattering rate cor- 
responding to a 10^^ cm~'^ donor concentration. Fig- 
ure 8 presents simulated electron velocity profiles in for- 
ward and reverse directions at the electric field profiles 
for 5 /im neck samples. For the sake of evidence, the 
electric field profiles taken from Fig. 6 are given, too. 
As one can see, the velocity overshoot/undershoot and 
strong nonlocality effects are clearly visible. 

The sharp peaks at the top of the velocity profiles are 
related to the locality in optical phonon emission rate 
Vop- Note that the velocity profile is asymmetric due to 
the presence of the different electric field gradients and 
exhibits nearly constant value in the forward direction 
ranging from 0.75 up to 3 /im. The Monte Carlo simula- 
tion for electric field profiles from Fig. 8 gives the coeffi- 
cient of current asymmetry of 0.48. The factor strongly 
depends on impurity scattering rate: it decreases as the 
scatering rate increases. The simulation with donor con- 
centrations 2x10^^, 1x10^^, and 5x10^^ cm""^ gives the 
coefficient values of 0.62, 0.48, and 0.25, respectively. 

To get close-up view in the physics behind, we have de- 
picted the velocity and Vop profiles in both, forward and 
reverse, directions (Fig. 9). Also, stationary velocity is 
added to illustrate the overshoot/undershoot areas. It is 
obvious that the electron velocity minima coincide with 
the Pop maxima forming thus "comb" in electron distri- 
bution. The effect is equivalent to the the free-carrier 
grating formation phenomenon caused by locality of op- 
tical phonon emission as already demonstrated in n^nri^ 
InN structuresJ^ 
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FIG. 9: (Color online) Monte Carlo calculations of drift veloc- 
ity and local scattering rate of electrons by optical phonons 
in "frozen" electric field in 5 /im neck sample. Note opposite 
phases in the velocity and the scattering rate. Background 
schematically shows the geometrical shape of the structure. 



V. POSSIBLE APPLICATIONS OF THE 
EFFECT 



Since the two-dimensional bigradient effect induces 
asymmetry in / — ^-characteristics, it is reasonable to 
consider its possible applications in detection of the high- 
frequency electromagnetic radiation. To test this idea, we 
have placed samples in a microwave field of 10 GHz fre- 
quency and recorded microwave- induced voltage arising 
over the ends of the sample. The results for two sam- 
ples of different neck width recorded at room and liquid 
nitrogen temperatures are displayed in Fig. 10(a). 

One can see that with the increase of incident power 
the signal increases linearly in all the range of studied 
power. Thus, due to the asymmetry in / — y curves 
caused by the two-dimensional bigradient effect the mi- 
crowave radiation can be detected^^. However, as es- 
timates show, the sensitivity at room temperature is 
rather low. For instance, it is only about 12 mV/W and 
40 mV/W for the samples with neck width of 12 /im 
and 2 /im, respectively. Since the effect is related to 
the electron heating, the decease in temperature should 
induce the increase in the detected signal. The tempera- 
ture effect, however, is not as large as might be expected 
from hot-electron physics where change of the signal is 
proportional to the carrier mobility. As one can see, at 
77 K, the detected voltage increases reaching the values 
of 40 mV/W (12 /im neck sample) and 150 mV/W (2 /im 
neck sample). In fact, the rise in only by a factor of about 
4, while the mobility increases about 14 times. We as- 
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FIG. 10: (Color online) Detected signal versus incident power 
at 10 GHz frequency in the samples fabricated from the struc- 
ture 2DEG-B. (a) - Data for two samples employing two- 
dimensional bigradient effect with neck sizes of 12 fim and 
2 ^m at room and liquid nitrogen temperatures, (b) - Com- 
parison of the two-dimensional bigradient effect-based sam- 
ples with bow-ties diodes containing 2DEG layer. Geomet- 
rical dimensions and are the same for both type of samples. 
Data recorded in 12 fim neck size samples at room and liq- 
uid nitrogen temperatures. Insets show measurement schemes 
and design differences. Note much higher sensitivity of the 
2DEG bow-ties diodes. 



sociate the outcome with a presence of bahistic electrons 
which cross the effective length of the samples without 
scattering Ji It is evident that for implementation of the 
effect into practical devices, such a sensitivity values are 
not sufficient. Following the ideology of the bigradient 
effect, one needs to produce much higher non-uniformity 
of the electric field in the sample's neck vicinity. This can 
be achieved easily by metallizing one of the asymmetric 
parts of the structure. This is evidenced by calculations 
of the distribution of the electric field given in Fig. 11. 

It is seen that the non-uniformity of the electric field 
with one part of the sample metallized is much stronger 
than in the conventional asymmetrically-shaped struc- 
ture. It is therefore reasonable to expect significantly 
higher values of the sensitivity. The experimental illus- 
tration is given in Fig. 10(b) where the bigradient struc- 
tures are compared with the so-called 2DEG bow-ties 
diodesl*— structures where one of the parts of the sam- 
ple is covered by metal. One can see that these devices 
are nearly in two orders of magnitude more sensitive com- 
pared to the structures of bigradient-design. 

Of particular interest for the direct applications are 
the values of the voltage sensitivity at room tempera- 
ture. In order to increase the sensitivity of the diode, we 
have chosen the design 2DEG-B (due to the higher value 
of the sensitivity at room temperature) and reduced the 
size of the apex down to 2 /im keeping the other geomet- 
rical dimensions the same (Fig. 12). As one can see, in 
this case the detected signals are significantly larger than 
for the 12 /im neck diodes at 300 K, and the sensitivity 
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FIG. 11: Calculated effect of the metallization of the B-part 
of the asymmetrically shaped sample on I — V characteris- 
tics. Room temperature value of the electron mobility fj,— 
4700 cmV(V-s) is used. The neck width is 5 /im. The filled 
squares denote the forward current //, the empty squares - 
the reverse current Ir- The coefficient of asymmetry is defined 

as {If-Ir)/If. 



amounts to 2.5 V/W, which is nearly one order of magni- 
tude higher compared to the diodes of 12 /zm neck width, 
where it is about 0.32 V/W. The decrease in tempera- 
ture down to 77 K allows one to increase the sensitivity 
up to nearly 40 V/W, i.e. close to a factor of 14.5 which 
corresponds to the change of 2DEG mobility. The latter 
effect - via thickness of spacer - is also nicely expressed 
in Fig. 12 were the detected signal in both devices, fab- 
ricated from 2DEG-A- and 2DEG-B-type structures is 
compared. 

Additionally, it deserves mentioning that the further 
decrease of the neck size down to 1 /xm - 800 nm, gives, 
however, no increase in the sensitivity (the data are not 
presented here). We attribute it to weaker coupling of 
the microwave field into the structure with the reduction 
of the neck width. 

As to dynamical range, the devices can be used to 
detect rather different power levels. As already in- 
dicated above, the bigradient modulation-doped struc- 
tures at 300 K in microwaves exhibit sensitivity of about 
12 mV/W and 40 mV/W for the samples with neck width 
of 12 /im and 2 /im, respectively, however, the sensing 
range is below 1 W in a pulse mode, i.e. the resis- 
tance to overloading is small. Bulk silicon-based 12 /im 
neck width structures are also not sensitive - only about 
0.1 V/W - but they can withstand up to 10 kW max- 
imum applied power at room temprature without any 
negative outcome. ^° Asymmetrically-shaped bulk GaAs 
structures containing n — n^-junction were found to be 
suitable for broadband sensing within GHz-THz frequen- 
cies at 300 K^ 

Finally, it is worth to remark on the operation of the bi- 
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FIG. 12; (Color online) Detected signal as a lunction ol inci- 
dent power at 10 GHz frequency microwave field in 2DEG 
bow-ties diodes of different geometry and design. Data 
recorded in 12 /^m and 2 /im neck size samples based on 
structures 2DEG-A and 2DEG-B at room and liquid nitro- 
gen temperatures. Linear function is also given as a guide for 
an eye. 



gradient structures within the THz range. One can note 
that at 300 K the sensitivity at 0.763 THz in the sample 
with neck width of 12 /im it is about 13 mV/W, i.e. close 
to the value in 10 GHz range. Narrowing the neck width 
down to 2 /xm, the sensitivity at 0.763 THz can increase 
up to 17 mV/W. Within frequencies 1.41 THz, 1.63 THz 
and 2.52 THz the detection observed is rather poor. We 
explain it by a weak couphng of the incident THz field 
into the structure. At hehuni temperatures, the sensitiv- 
ity reaches the value of 60 mV/W at 0.763 THz for the 
12 /Ltm neck width sample.^— 
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FIG. 13: Change of the relative current in 5 /im neck size 
sample at liquid helium and room temperatures. Insert de- 
picts the same plot for comparison in n— doped asymetrically- 
shaped silicon samples with specific resistance of 2.5 f2-cm at 
300 K. The length of the sample is 23 mm, angles of the asym- 
metry are ai=45'' and 02 = 13°. More details can be found in 
Ref. 19. 



in 5 /im neck size sample at 4.2 K and 300 K. For com- 
parison of the features, inset shows equivalent bulk data 
measured in n-type silicon according to Ref. 19. 

It is seen that in the bulk semiconductor the nonlocal- 
ity is associated with the sign change in the asymmetry 
clearly observed within 1-10 V voltage scale at 77 K (see 
inset). In contrast, in modulation-doped structure, no 
such peculiarity is observed. It is worth noting, however, 
that above 0.3 V the change of the relative current is 
nearly independent of voltage - as was indicated earlier 
(Fig. 7, above 0.25 V range, where asymmetry of / — T^ 
characteristics strongly increases) - we attribute it to the 
apperance of the nonlocal drift velocity. 



VI. TWO-DIMENSIONAL VS. BULK 
BIGRADIENT EFFECT 

The effect was originally discovered by measuring 
I — V characteristics in bulk semiconductors, n-Ge and 
n-Si, at room and liquid nitrogen temperatures .-SJ^ In 
what follows we compare the two-dimensional bigradi- 
ent effect observed in asymmetrically in-plane shaped 
GaAs/AlGaAs modulation-doped structures with its 
bulk equivalent desribed in details. ^^ 

The inherent fearure of both - the two-dimentional and 
bulk - effects is asymmetry in I — V characteristics. Dif- 
ferent material parameters and geometrical dimensions 
of the samples determine rather distinctive voltage scale 
and temperature conditions for experimental observa- 
tion. The effect of the dimensionallity we associate with 
unusual manifestation of the nonlocality in drift veloc- 
ity. Figure 13 presents a change of the relative current 



VII. CONCLUSIONS 

Asymmetrically in-plane shaped GaAs/AlGaAs 
modulation-doped structures of various design are stud- 
ied within 4-300 K temperature range. We have shown 
that current -voltage characteristics of such structures at 
low, 4-80 K, temperatures exhibits pronounced asym- 
metry which we explain by a two-dimensional bigradient 
effect. This phenomenon can be decomposed into two 
constituents - carrier accumulation and nonlocality in 
the electron drift velocity - each of which manifests 
itself depending on the values of the in-plane electric 
fields and their gradients. Varying these parameters 
via change of the sample neck we have inferred that 
the effect is induced by different accumulation of 
two-dimensional electrons due to asymmetrical shape 
of the structure when the neck size down is 12 /im; 
in 5 /im neck size structures in electric fields close to 
3 kV/cm, the nonlocality of the electron drift velocity 



becomes predominating causing thus further increase 
of the asymmetry in / — T^ curves. The theoretical 
models suggested conform well with the obtained ex- 
perimental results. Monte Carlo simulation has shown 
strong nonlocal effects in the electron drift velocity and 
presence of overshoot/undershoot regions depending on 
the electric field and its gradient values. The features 
of the two-dimensional bigradient effect and its bulk 
equivalent are compared and possible applications of 
this phenomenon to detect electromagnetic radiation of 
GHz and THz frequencies are discussed as well. 
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